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MULTI-ZONE  MODELING  OF  IMPURITY  REDISTKIHUTION 


IN  ION-IMPLANTED  MATERIALS 

by 

Richard  Kwor 

ABSTRACT 

Implanted  impurity  redistribution  has  been  observed  during  annealing  of 
many  ion- implanted  materials.  Experimental  evidence  suggests  some  position 
dependence  in  the  redistribution  process.  The  tail  region  of  ion-implanted 
impurity  profiles  usually  exhibits  faster  diffusion  than  the  near-surface  re¬ 
gion.  In  this  paper,  a  multi-zone  model  for  the  redistribution  of  implanted 
impur it ies  is  presented .  The  implanted  substrate  is  considered  as  a  stra¬ 
tified  medium  with  zones  where  a  local  diffusion  equation  is  obeyed,  and  an 
effective  diffusion  coefficient  is  defined  within  each  zone.  The  basic  formu¬ 
lation  of  the  model  and  its  mathematical  background  are  discussed.  The  multi¬ 
zone  equations  are  solved  using  the  Crank-Nicol son  method.  A  computer  program 
is  used  to  generate  a  plot  of  the  post-annealing  redistributed  impurity  pro¬ 
file.  The  model  is  applied  to  the  case  of  sulfur-implanted  GaAs  for  dose 
range  of  4  x  10^^  to  4  x  10^^  cm**^^  with  energies  of  120  keV  and  300  keV. 

Good  agreement  is  obtained  between  the  computer  generated  profile  and  the  SIMS 
experimental  profile. 
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INTRODUCTION 


Ion  implanLat Ion  offers  a  number  of  advantages  in  the  fabrication  of 
semicijnduc tor  devices  and  is  now  an  established  technology  for  production. 

Its  attractiveness  comes  from  the  ability  to  achieve  doping  profiles  which  are 
controllable,  reproducible  and  sometimes  not  easily  obtained  by  other  tech¬ 
niques,  In  the  designing  and  characterization  of  ion  implanted  devices,  it  is 
therefore  very  important  to  know  or  to  be  able  to  predict  the  doping  profile. 
In  most  cases,  the  final  carrier  concentration  profile  is  closely  related  to 
the  implanted  impurity  distribution  and  the  post -imp Ian  tat  ion  anneal  ing  . 
High-temperature  annealing  is  necessary  for  the  removal  of  radiation  damages 
and  the  activation  of  implanted  atoms.  However,  the  diffusion  of  the  impurity 
atoms  during  annealing  often  results  in  substantial  movement  of  these  atoms, 
altering  the  impur ity  distribution.  The  pre-annealing  impurity  distribution 
is  usually  described  by  the  profiles  of  LSS  (1),  joined  hal f-Gaussian  (2),  and 
more  recently*  the  Pearson  IV  (3,  A),  DepiMiding  on  the  implantat ion  condi¬ 
tions,  these  profiles  have  various  degrees  of  success.  For  post-annealing 
prof  i  les  ,  however ,  the  ef  feet  of  d i f fus ion  must  general  ly  be  inc  luded  .  The 
commonly  used  model  for  impurity  diffusion  during  annealing  assumes  an  initial 
Gaussian  distribution  and  a  constant  diffusion  coefficient  (5,  6).  In  the 
case  of  capped  substrate  such  as  Si3N4  on  GaAs ,  the  encapsulant  is  sometimes 
considered  a  continuation  of  the  substrate.  According  to  this  model,  anneal¬ 
ing  will  result  in  a  broader  and  lower  Gaussian  distribution.  This  type  of 
diffusion  can  only  be  regarded  as  an  approximation,  whereas  the  actual  diffu¬ 
sion  of  implanted  species  is,  in  general,  far  more  complex.  Factors  such  as 
lattice  defects  may  play  an  important  role  in  the  diffusion  mechanics  of  im¬ 
purity  atoms.  During  implantation  different  types  of  defect  occur  in 


damage  clusters  around  the  path  of  the  implanted  particle  (7),  The  nature  and 
spatial  distribution  of  radiation  damage  is  I'xtreinely  complicated.  'Ine  exact 
class  of  defect  or  defect  complex  depends  on  many  factors,  including  the  mass 
and  energy  of  the  inc  ident  part ic les  and  the  mass  of  the  target  atom.  Gen¬ 
erally,  the  damage  profile  does  not  coincide  with  the  implanted  ion  profile; 
the  degree  of  disparity  between  these  two  being  a  function  of  the  mass  differ¬ 
ence  between  the  bombarding  ion  and  the  lattice  atoms.  Moreover,  when  the 
solid  solubility  of  the  implanted  spec ies  in  the  substrate  is  exceeded ,  as  in 
the  case  of  high  dose  implantation,  precipitates  sometimes  occur,  thus  chang¬ 
ing  the  local  characteristics  of  the  medium.  Clearly,  the  pre-annealing  sub¬ 
strate  cannot  be  considered  as  a  homogeneous  medium  for  diffusion  of  the  im¬ 
planted  impurity  atoms.  Another  important  parameter  which  may  also  have  a 
considerable  influence  on  the  final  distribution  of  implanted  ions  is  the 
subst rate-encapsu I  ant  boundary  condition.  Most  encapsulants  behave  as  bar¬ 
riers  to  outd  i  f f us  ion  ,  leading  to  a  build  up  of  impurity  atoms  at  the  sub¬ 
st  rat  **  surface,  but  it  is  also  possible  ft>r  some  to  act  as  sinks.  Tlie  strain 
at  the  substrate-encapsulant  interface  also  has  an  effect  on  the  diffusion  of 
the  implanted  imparities  (8).  This  complexity  is  supported  by  the  experimen¬ 
tal  evidence  that  many  post-annea I  ing  impurity  profiles  cannot  be  accurately 
described  by  the  simple  diffusion  theory.  Examples  of  long  diffusion  tails 
are  numerous  (9,  10,  ll).  They  suggest  that  the  implanted  medium  may  be  char¬ 
acterized  by  two  or  possibly  more  zones  with  different  diffusion  constants* 
This  results  in  a  multi-zone  modeling  scheme  for  the  ion  implanted  medium. 

The  medium  is  divided  into  zones  where  a  local  diffusion  equation  is  obeyed 
with  an  effective  diffusion  coefficient  which  lumps  the  effects  of  the  defects 
and  possible  prec ipi tates .  With  appropriate  boundary  and  interface  condi¬ 
tions,  this  model  can  be  used  to  more  accurately  describe  the  post-annealing 


impurity  profiles.  In  its  extended  form  it  ean  be  used  to  model  out d i f fus i on 


of  implanted  species  into  the  encapsulanl.  In  this  paper,  the  basic  formtila- 
t  ion  of  the  model  and  its  mathematical  background  are  iliscussed.  A  computer 
program  is  then  used  to  generate  a  plot  of  the  post-annealing  redistributed 
impurity  profile.  Reasonable  agreement  is  obtained  between  the  computer  gen¬ 
erated  profile  and  the  SIMS  atomic  profile  of  sol  fur- implanted  Cr-doped  GaAs 
for  dose  range  of  4  x  10^^  to  4  x  loO 

cm  ^  with  energies  of  120  keV  and  300 

keV.  This  provides  a  simple  and  yet  more  accurate  way  of  estimating  the  dif¬ 
fusion  coefficients  of  an  ion  implanted  medium.  The  advanced  matrix  formula¬ 
tion  of  the  model  allows  the  computer  program  to  be  independent  of  the  initial 
profile  and  also  virtually  independent  of  constrictions  in  the  boundary  condi¬ 
tions.  The  method  is  thus  very  suitable  for  profile  analysis  in  ion  implan t a- 
t ion  when  the  initial  profiles  cannot  be  defined  by  analytical  forms  as  in  the 
case  of  non -Gauss ian  profiles  or  mu Itiple  implantations. 
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BASIC  FORMULATION  OK  MULTl-ZONK  MODEL 


The  basic  concepts  of  a  double- zone  m<idel  will  first  be  discussed.  In 
later  sections,  these  basic  concepts  are  ^generalized  for  the  multi-zone  case. 
P’ig.  I  shows  the  geometry  of  a  double-zone  model.  The  general  continuity 
equation  for  the  position  dependent  diffusion  cot^fficient  is: 


3C 

3t 


3 

3x 


3  c 

D(x)  - 

3  X 


[11 


where  C(x,  t)  is  the  concentration  of  the  implanted  impurity  atoms,  D(x)  is 
the  position  dependent  diffusion  coefficient,  t  is  the  time  and  x  is  the 
depth.  For  a  double-zone  model,  each  of  the  two  zones  is  assumed  to  have  a 
constant  diffusion  coefficient: 


D(x) 


for  0  ^  X  ^  a 
for  a'*’  <  X  < 


(2) 


where  D|  and  D2  are  the  effective  diffusion  coefficients  for  zone  I  and  zone  2 
respectively,  a  is  the  location  of  the  zone  interface  and  ^  is  the  location 
of  the  second  zone  boundary.  Thus,  Eq .  [1]  becomes  the  familiar  Pick’s  second 
I  aw  inside  each  zone : 
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for  0  <  X  <  a  (3a] 


for  a  <  X  <  it  [3bl 


where  C\  and  C2  are  the  concentrations  in  their  respective  zones.  If  the 
initial  condition,  i.e.  the  as-implanted  impurity  profile  is  F(x),  then 
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for  0  <  X  <  a  (4al 


for  a  <  X  <  i  (^b) 


^ i  — "  +Pi(x)Ci=fi(x,t)  [5] 

3  X 


at  the  boundary  of  the  i th  zone  where  the  concent  rat  ion  C{(x,  t)  ic  defined. 


The  function  fi(x,  t)  is  seen  as  a  source  of  impurities  located  at  the 
boundary.  The  function  P|(x)  is  usually  a  constant  related  to  the  absorption 
or  segregation  properties  of  the  boundary.  If  p£  =  0,  then  Eq.  (5]  reduces 


to 


3C; 

-  =  f i(x,  t)  (6) 

3  X 


For  a  reflective  boundary,  i.e.  no  diffusing  atom  crosses  it,  =  0,  and 

3  X 

there  is  no  net  flux  across  the  boundary.  And  if  q|  =  0  and  p^  ^  0,  then 


PiC2®=fi(x,t)  [7] 


Ci(x)  =  F(x) 
C2(x)  =  F(x) 

The  general  boundary  condition  is 

3C; 


An  example  for  this  case  is  C£(x,  t)  =  C}^,  the  background  concentration  of 
impi»rity  in  the  substrate  when  I  is  very  large. 

In  our  simple  double-zone  model,  three  boundary  conditions  will  be 
assumed.  A  perfectly  reflective  cap  is  assumed  to  exist  at  the  surface. 
Thus , 


3Ci 

3x 


=  0 


[8] 


I  X  =  0 

For  large  i  ,  it  is  reasonable  to  assume  that  no  diffusing  impurity  will  cross 
the  boundary  at  x  *  &  for  the  short  annealing  time  used  In  ion  implantation. 
Then 


3 


3  X 


X  =  ^ 


0 


19] 


And  at  tho  zoiio  iaL».»rtaco  of  x  =  a,  the  flux  in  /ione  1  t*  qua  Is  that  in  zone  2. 
Hence  , 


n 


I 


dC\ 
3  X 


3  X 


at  X  =  a 


Iiol 


Physically,  this  means  that  the  inter  zone  boundary  is  nonabsorbing  and  dot's 
not  segregate  impurities.  Also,  in  this  model,  the  locations  of  the  ziine 
boundaries  are  assumed  to  be  t  ime- independen  t  .  'then ,  given  an  initial 
Condition  C(x,  0)  =*  F(x),  and  with  the  boundary  conditions  specified  by  the 
Fqs.  [^],  (9),  [10),  Eq .  [3]  can  be  solved  to  yield  the  functions  C\  and  C2 
for  a  specific  combination  of  annealing  time  and  the  effective  diffusion 
C4)e f f  ic  ien t s  in  the  two  zones.  The  mathematical  method  is  outlined  in  the 
next  sect  ion , 
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DISCRETI/ATION  OF  THK  DIFFUSION  F()nATl()NS  FOR  DOUBLF'-ZONE  MODEL 


To  K(|s  .  1  la  1  aiul  |  D'*  ]  lor  tho  Uoub  I  oast*,  the  eonl  iiiuons  set 

of  space- time  points  {(x,  t)}  pairs  is  substituted  with  a  of  discrete 

points  {  =  nA  X  ;  t,-j  =  nA  t}  ,  where  Ax  and  At  are  small  intervals  of  distane 

anil  time  respectively,  m  and  n  are  non-negative  integers,  in  =  0 ,  1,  2,  ...»  L; 

n  =  0,  1,  2,  T.  At  the  beginning  of  first  zone,  x^,  =  0,  m  =  0;  and  at 

the  end  of  second  zone,  ,  m  =  L.  The  total  annealing  time  is  equal  to 

TA  t ,  The  func  t  ion  C(x,  t)  at  the  grid  pi^'int  (m,  n)  is  approximately 

n 

represented  by  the  point  Cjn,  Using  the  Crank-Nico Ison  method  (12),  Eq .  [3] 


c an  be  transformed  to : 


fii  .  pn  + 1 

,.  n  +  ^-  ^  -rg  > 


p-'  C"  =■  Di  O  +  p  -  (  - 

t  m  XX  2 


where  is  the  diffusion  coefficient  of  the  ith  zone. 


0  -  the  finite  difference  operator  for  the  forward  time  step, 

t  f.n+1  n 

p  >  cn  = 
t  m  At 


p  ^  ~  the  finite  difference  operator  for  the  fcirward  distance  step. 

^  — - 

X  m  Ax 


the  finite  difference  operator  for  the  backward  distance  step. 

-  cn  =  SeJLSezL 

X  m  Ax 


Eq .  (Ill  can  be  expanded  to  a  set  of  equations  (13): 


n+ 1  A  :  n+ 1  n+ 1  n  A  ;  n  n 

(l+Xi)C  -  -  (C  +C  )  =  (l-Xi)C  + -  (c  +c  ) 

m  2  1  ni*  i  m2  ^  ^ 
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L 

wher<!*  ^  2  ,  thi*  normal  <iiffusiof^  ooeffioiont  for  tlie  il.h  7.oti»*. 

(Axl? 

Tht?  initial  contliiion  (n=0)  is  roadily  iiu'or poral4*d  in  this  system,  and  Lht* 
boundary  conditions  in  the  discrett*  space-grid  art*; 

(1)  At  the  beginning  of  the  first  /.one,  Xq)  “  ^  ^  •  !12]  thus 

becomes 


( I  >  X i)C, 


n+1  n+i  ti  *  n  n 

(C  +  C  )  =  (1  -  X  1)0^3  ^  ±A  iC  +  C  ) 


If  we  assume  a  perfectly  reflecting  surface,  then 


lx  =  0 


which  in  the  discretized  format  becc^mes 


o  n 

0  c 

X  m 


^ m+ 1  ^m“l 


there  fore 


iP  t  + 


n  n 
Cl  -  c_i 


n  n 

or  C  =  C 
I  -1 


and  Kq .  (13 1  bec«)nies 


n+I  n+I  n  n 

(i+xpc  -XiC  ^  [\  -  \\)  C  +XiC 
o  I  o  I 
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(2)  At  the  end  of  the  second  zone,  in  =  L,  i  =  2, 

o  n 

l)  C  =  0,  Eq ,  [  1 2 1  becomes 

X  o 


n+1  n+l  n  n 

(1+X2)C  =(1“X2)C  +^2^^ 

L  L-1  L  b-1 


(3)  At  the  zone  interface,  =  n,  m  =  A,  and 


[15] 


\)[  0 

X 


n 

C 

A 


then  the  inter-zone  equation  can  be  shown  to  be; 


X  1  \  2  1  X  1  n+ 1  X  9  n+1 

(1  +  +  )  C  C  C 

2  2^2  2 


Xi  \9  n  Xi  n  Xon 

(I  -  —  )  c  +—  c  c 

2  2  ^2  A-I2  1 


116) 


For  m  =  0,  1,  2, . ,L,  the  system  of  equations  Eq .  [12]  can  be  represented 

by  a  matrix  equation 


A  C"'*’!  =  B  C’ 


[17] 


where  = 


•  • 

r  *1 

n+1 
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C 

and  = 

C 

0 

0 

n+ 1 

n 

C 
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DepentlLiiR  on  the  surface  contJition  (at  =  0),  the  in  ter  zone  conditions  and 

the  last  zone  boundary  condition  (at  x,j,  -  i),  the  A  and  U  matrices  will  take 

slightly  different  forms,  A  computer  program  is  then  used  to  solve  the  matrix 

n 

e(|uai  ion  Fq ,  (I7l  and  calculate  the  valuos  of  C  at  difterent  grid  points 

m 

(m,  n).  The  red  i  st  r  ibuted  pr<>file  after  a  certain  annealing  time  can  then  be 
plotted. 
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THK  MULTl-ZONfc:  MODHL  WITH  OiNKKAl.  BOUNDARY  CONDITIONS 


The  multi'zotu*  mi)de  I  with  bounJ.iry  coiiii  i  t  i  on  s  (  Eq  .  (5))  follows 

the  same  arguments  as  the  double-zone  model.  Inside  the  ith  zone,  where  i  = 
1,  2,  1, . ,1,  Kick’s  law  is  assumed  to  be  i)l)eyed : 


-  =  Di  - 

3t  3^2 


ins  ide  ith  zone 


Assume  the  b<>undary  condition  at  the  surfact?  to  be 


+  PjCi  =  4>  1  (t) 


at  X  =  0 


And  at  X  =  It  the  boundary  condition  is: 


+  P I C  [  =  |}>  I  ( t ) 


Her.'  1  ( t. )  aiul  are  time  dependent  sources  (or  sinks)  at  the  first  and 

last  boundaries.  The  functions  P\  and  Pj  art*  relatt*d  to  the  permeability  of 
the  bt^undaries.  These  boundary  conditicjns  ar»*  ru>w  discretized: 


os  s  s 

n  C  +  P|C  =  <p 
X  o  o  I 


os  s  s 

D  C  +  p  C  =  <j) 

X  L  T  L  I 


Xj^  =  i 


where  s  =  n ,  n+1 


The  first  and  last  equations  in  Eq,  [12]  then  become : 


n+l  n+1  n+1 

(l  +  (l-Axp|)A]]C  -A|C  +A|Ax^ 
o  1  I 


n  n  n 

[1  -  (1  -  Axpj)AilC  ♦A^C  -Ai 

o  1  I 


[251 


I  n  I  11  f  I 

[  I  +  (  1  +  A  xpT  )A  i]C  -AtC  -AtAx'> 

L  f,-l  I 


n  (1  II 

=  1  I  -  (  I  ^  A  xdt  )A  r )  C  ♦  a  t  c:  ♦  a  t  a  X({> 

L  t,-!  I 


[26) 


Eqs.  (25)  and  (26)  are  seen  to  reduce  to  Eqs .  (14)  and  (15)  for  the 
double-zone  case  i f  1  *  2 ,  Pj  *  Pj  =  0  and  ^  f  >0,  For  the 

multi-zone  case,  another  matrix  equation  can  be  obtained: 

A  =  B  [27! 


s  s  T 

where  ^  ^  =  (A  |A  ^  0  0 . A  |A  x4>  ]  s  -  n ,  n  +  1 

If  the  inter-zone  conditions  are  kept  the  sanu^  as  that  for  the  double-zone 
cast%  i.e.  assuming  no  interface  resistance  or  enhancement  to  particle  flow, 
approaches  similar  to  that  for  the  double-zone  model  can  be  used  to  find  the 
matrices  A  and  H.  And  again  the  red  i  st  r  i  hut  t*d  profiles  can  be  calculated 
given  the  annealing  time  and  the  set  of  diffusion  coefficients  for  the  I 

z.Jfies  . 
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APPLICATION  OK  THK  MODKl.  TO  THK  KKI)  I  S  TK I KUTION  OK  SULKUK 
IN  ION-IMPLAN7>:i)  (LiAs 

In  this  st»c  t  ion  the  inu  It  i  -  zone  moile  1  dese  r  i  bed  in  the  prev  i  ons  s<*c  t  i  ons 
is  applied  to  some  pri)files  sulfur  implanted  in  <100^  oriented  Cr-doped 
sem  i  -  i  n  so  I  a  t  i  ng  OaAs  substrate,  Th»*  iinpl  antat  ions  were  done  at  room  tempera¬ 
ture  and  in  a  non-channe 1  ing  direction  to  doses  ranging  from  U  x  10*^  to 

4  X  loO  cm''^ .  After  implantation,  the  samples  were  capped  with  either  Si3N4 
or  Si02  and  then  annealed.  There  is  a  significant  redistribution  of  sulfur  as 
a  result  of  diffusion  during  annealing.  As  indicated  by  the  SIMS  prcjfiles 
(Figs,  2-5),  the  as- implanted  profile  cannot  he  be  adequately  described  using 
t,SS.  Pearson  IV  distribution  was  therefore  us**d  as  the  initial  condition  for 
the  multi-zone  model.  Using  the  same  appr^jach  as  that  in  Pef,  (A),  the  four 
•nomeits  for  the  Pear  soil  IV  distribution  were  ftujiul.  fhesi*  four  tnrjments  were 
then  fed  to  the  computer  which  generates  the  initial  condition  for  the  model. 
Attempts  to  match  the  pos t -annea I i ng  profiles  for  these  sulfur  implants  with 
Pears.jn  IV  distribution  have  been  unsuccessful  because  of  the  long  diffusion 
tails  involved.  However,  as  will  be  seen  later,  the  triple-zone  model  can  be 
us»?d  to  adequately  describe  these  profiles.  Double-zone  model  has  also  been 
tried,  hut  it  does  not  give  as  good  results.  As  discussed  earlier,  the  rmulel 
accepts  initial  condition  of  any  distr  ibut  i»)n .  Pearson  IV  has  been  used 
throughout  because  it  gives  satisfactory  results  for  these  particular  appIica- 
t  ions.  From  the  SIMS  data  shown  in  Kigs.  2-5  the  different  zones  in  the  re¬ 
distributed  profiles  are  apparent •  Depend ing  on  the  anneal  ing  t ime ,  tempera¬ 
ture,  implantation  energy  and  dose,  the  effective  diffusion  coefficients  in 
different  zones  are  different.  In  Figs.  2-4,  the  impurity  concentrations  near 
the  surface  exceed  the  solid  solubility  due  to  high  dose  implantation.  Pre¬ 
cipitation  and/or  a  large  amount  of  unannealed  ilisorder  in  this  region  results 
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in  n  vt*ry  small  diffusion  coe  f  f  i  c  i  l  .  l)t*i»per  in  the  substrate,  the 

s  til  fur  atoms  diffuse  with  a  larp,er  effective  diffusion  coefficient.  Under 
such  circumstances,  significant  redistribution  of  sulfur  only  occurs  In  the 
tail  region  of  the  profile.  Tn  lower  dose  implants,  there  is  a  much  smaller 
<!ifference  in  the  effective  diffusion  ctie  f  f  ic  lent  s  of  different  zones,  except 
in  the  lower  temperature  annealing.  Fig,  S  gives  such  an  t*xample. 

The  multi-zone  model  was  applied  to  four  cases.  In  each  case,  the 
Pearson  IV  moments  for  the  as-implanted  profile  were  found  from  the  SIMS  data 
and  used  to  generate  the  initial  profile.  The  redistributed  profile  after 
anneal  ing  was  then  calculated  using  a  computer.  The  four  cas<^s  are  discussed 
be  1 ow , 

The  SIMS  unannealed  and  annt»aled  profiles  of  Cr-doped  GaAs  implanted  with 
1  20  keV,  ini'*  cm’  -2  32s  are  shown  in  Fig.  2.  These  profiles  were  reported 
earlier  in  Ref ,  9.  The  annealed  profile  shows  a  very  distinct  mult i-zone 
character  ist  ics .  Pearson  IV  disLr ihut  Lon  was  used  to  describe  the 
as  -  imp  I  aji  t  eil  profile  which  is  seen  tt»  deviate  substantially  from  bSS.  With 
the  choice  of  ft)ur  moments  of  Rp  =  O.l  pm,  o  ==  0.08  prn,  Y  “  “1.7  and  3  =  AO,  a 
vood  fit  was  obtained  except  near  the  tail  region  whert*  the  Pearson  IV 
iistribution  has  a  slightly  sliallower  penetration.  After  a  IS  minute  anneal 
at  900“C,  performed  with  an  Si^N/^  encapsulant,  a  pri^nounced  tail  appeared  in 
the  annealed  profile.  The  tail  is  seen  to  start  at  about  2  x  cm*"^,  which 

is  ci>mm»>nly  accepted  as  the  solid  solubility  of  implanted  sulfur  in  GaAs,  To 
model  this  annealed  profile,  a  triple-zone  scheme  was  used.  The  Pearson  IV 
Iistribution  with  the  above  mentioned  four  moments  was  used  as  the  initial 
•riindition.  The  triple-zone  profile  with  interfaces  at  0.2  Mm  and  0.22  pm,  and 
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three  effect  ive  diffusion  ct)e  f  f  i  c  lent  s  of  Uj  -  4  x  lO”^^  cin^/sec,  D2  ~  8  x 
10"®^  om^/  sec ,  D3  =  1/56  x  10”^^  cm^/  sec ,  is  shown  in  K  •  6 . 

In  Fig.  3,  the  SIMS  unannenled  and  annealed  profiles  of  Cr-doped  GaAs 
implanted  with  120  keV,  10^^  cm”^  are  shown.  The  f(jur  fn<jments  used  in  the 

Pearson  IV  distribution  for  the  as-implanted  profile  were  Rp  =  O.ll  Urn,  0  = 
0.097  lim,  Y  =  *2.5  and  3  =  40.  Again  p  a  good  fit  was  obtained  except  near  the 
tail  region  where  the  Pearson  IV  has  a  shallower  penetration,  which  is  common 
in  all  the  cases  reported  in  this  work.  After  15  minutes  anneal  at  900*C, 
performed  with  an  Si3N4  encapsulant ,  a  pronounced  tail  similar  to  the 
above-discussed  10^^  cm”^  dose  case  appeared.  The  triple-zone  profile  for  the 
annealed  curve  is  shown  in  Fig.  7.  The  interfaces  used  were  at  0.21  pm  and 
0.23  pm.  The  three  effective  diffusion  coefficients  were  Dj  =*  4  x  10”^^ 
cm- /sec,  D2  =  8  x  10“^^  cm^/sec  and  D3  =  3.6  x  10”^^  cm^/sec. 

The  SIMS  data  studied  in  the  last  case  here  were  also  from  Ref.  10. 
Several  annealed  curves  for  different  temperatures  are  shown  in  Fig.  5. 
Multi-zone  character  ist  ics  of  the  profile  are  more  apparent  in  the  lower 
temperature  annealed  curves.  Triple-zone  modeling  was  performed  fur  the  700*C 
annealed  prcjfile.  The  SIMS  as-implanted  profile  from  Fig.  5  was  first 
calihraled  in  terms  of  concent  r  at  ion  by  integrating  the  area  under  the  curve 
and  setting  it  equal  to  the  area  under  the  700‘’C  annealed  curve.  The  four 
moments  of  the  Pearson  IV  were  then  found:  Rp  =  0.24  pm,  0  =0.12  pm,  Y  =  “*1 
and  3  =  40.  This  initial  condition  was  then  used  for  the  triple-zone  model. 
(Fig.  9).  The  interfaces  used  were  at  0.42  pm  and  0.60  Pm.  The  effective 
diffusion  coefficients  were  I)|  =  9  x  10"^^  cm^/sec,  D2  =  2.07  x  10*^^  cm^/sec 
and  D3  =  2.7  x  10“^^  cm^/sec. 

The  above  results  indicate  that  for  high  dose  implants,  there  is  very 
little  diffusion  in  the  first  two  zones,  while  the  effective  diffusion 
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CDefficient  in  the  thir*.i  zot\e  is  boCweeii  10*^^  .ind  10“^^  cm^/sec,  the 

valiies  of  sulfur  diffusion  coefficient  r^'porled  by  otfier  researchers  (6,  14, 
15,  16).  For  lower  dose  and  hi^'h  energy  implants,  the  difference  in  diffusion 
coefficients  of  differei\t  7,ones  depends  on  the  annealing  temperature.  At 
annealing  temperatures  below  700°C,  there  is  little  diffusion  in  the  first 
but  the  diffusion  in  the  tliird  zone  is  quite  fast. 


SUMMARY  AND  DISCUSSION 


In  all  casi»s  sliulioci,  reasonable  fit  was  obtained  between  the  triple-zone 
profile  and  the  cor respond ing  SIMS  annealed  curve,  with  best  fit  obtained  in 
the  first  and  third  zones.  Since  there  is  a  change  of  as  much  as  several 
orders  of  magnitude  in  the  effective  diffusion  coefficients  from  the  first  to 
the  third  zones ,  the  second  zone  is  in  reality  a  transition  zone  in  which 
there  is  a  continuous  change  of  diffusion  coefficient.  Because  of  the 
piece-wise  linear  nature  of  the  triple-zone  model,  only  one  effective 
diffusion  coefficient  is  chosen  for  the  second  zone.  Thus  the  greater  the 
difference  between  the  diffusion  coefficients  in  the  first  and  third  zones, 
the  m<>re  difficult  it  is  to  have  an  accurate  fit  of  the  profile  in  the  second 
zone.  The  results  from  the  multi-zone  modeling  show  tliat  the  impurity 
diffusion  is  related  to  the  implantation  energy,  dost*  and  annealing 
temperature.  It  is  believed  that  at  different  implantation  energy  and  dose, 
different  damages  are  created,  leading  to  different  effective  diffusion 
coefficients  in  the  third  zone.  Higher  implantation  energy  results  in  heavier 
damage  and  a  larger  effective  diffusion  coefficient  in  the  third  zone.  Larger 
doses  have  a  similar  effect.  On  the  other  hand,  the  heavy  damage  near  the 
surface  from  higher  implantation  energy  or  larger  dose  results  in  a  slow 
diffusion  in  that  region.  Fig.  5  gives  a  good  indication  of  the  influence  of 
annealing  temperature  on  the  difference  in  effective  diffusion  coefficients  of 
different  zones.  In  this  case,  when  the  implantation  energy  is  high  and 
antiealing  temperature  is  low  (around  700*C),  not  much  the  implantation 
damages  near  the  surface  are  removed  during  annealing,  and  there  is  a  great 
difference  in  the  effective  diffusion  coefficients  of  the  different  zones.  At 
h  igher  anneal  ing  temperatures ,  the  implantat ion  damages  are  removed ,  the 
difference  in  diffusion  coefficients  in  different  zones  is  reduced.  At  850*C 


tempo ra t  iir^^ ,  the  annealed  profilt*  has  almost  lost  its  multi-zone 
characteristics. 

In  summary,  the  multi-zone  model  ^ives  a  rt' a  son  ably  accurate  description 
of  the  post-anneal  inp  atomic  profile  f*)r  ion- impl  an  tt^d  impurities.  It  also 
j>ives  some  insight  of  the  physical  process  occuring  during  the  post- 
implantation  annealing,  yielding  some  in  format  it)n  of  the  position  dependence 
of  the  diffusion.  Using  the  model,  effective  diffusion  coefficients  in  d i f- 
ferent  zones  of  the  implanted  material  can  be  mure  accurately  calculated. 

This  model  accepts  any  profiles  as  initial  condition,  and  is  thus  very  suit¬ 
able  for  profile  analysis  in  ion  implantation  when  the  initial  profiles  cannot 
he  defined  by  analytical  forms  as  in  the  case  of  non-Gauss i an  profiles. 
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KKUiKR  CAPTIONS 


Fig.  1.  Configuration  of  a  tlouhle-zono  model  for  impurity  redistribution. 

Fig.  2 .  SIMS  unannea  I  ed  and  atinea  1  ed  profiles  of  Cr -doped  GaAs  impl  an  t»*d  with 

120  keV,  10 15  mimite  .iniu>al  was  performed  with  an 


Fig.  3. 


Fig.  4. 


Fig.  3. 


Fig,  6 . 


Fig.  7. 


Fig.  8. 


Si3N4  enc a psul ant. 

SIMS  unannealed  and  annealed  profiles  of  Cr-doped  GaAs  implanted  with 
120  keV,  10^5  cm*  ■2  32s.  15  minute  anneal  was  performed  with  an 

Si3N4  encapsulant. 

SIMS  unannealed  and  annealed  profiles  of  Cr-doped  GaAs  implanted  with 
300  keV,  4  x  10^^  cm“2  ^^S.  20  minute  anneal  was  performed  with  an 

Si02  encapsulant •  (After  C.  A,  Evans  et  al . ,  Ref.  9) . 

SIMS  unannealed  and  annealed  profiles  of  Cr-doped  GaAs  implanted  with 
300  keVj  4  X  10^^  cm“^  ^^S.  20  minute  anneal  was  performed  with  an 

Si02  encapsulant ,  (After  C.  A.  Evans  et  al . ,  Ref.  9) • 

SIMS  profile  of  120  keV,  IQi^  cm“2  implant  annealed  at  900°C  ff»r 

15  minutes;  Pearson  IV  profile  initial  condition;  and  triple-zone 
profilf^  with  interfaces  at  0.2  vim.  0,22  v^m  and  02,  O3  ennal 

4  X  8  X  1.36  X  cm^/sec  respectively. 

SIMS  profile  of  120  keV,  10^^  cm"^  ^^S  implant  annealed  at  900°C  for 
15  minutes;  Pearson  IV  profile  initial  condition;  and  triple-zone 
profile  with  interfaces  at  0.21  pm,  0.23  pm  and  Dj,  D2,  D3  equal 
4  X  8  X  10“^^,  5.6  X  10“^^  cm^/ sec  respect  i vely . 

SIMS  profile  of  300  keV,  4  x  10^^  cm*"^  implant  annealed  at 

840®C  for  20  minutes;  Pearson  IV  profile  initial  condition;  and 
triple-zone  profile  with  interfaces  at  0.39  pm,  0.60  pm  and  Dj,  D2, 

D3  equal  9  x  1 0*"  ^  ^ ,  1.44  x  1  ^ ,  1.98  x  1 0"*^  ^  cm^/  sec  re  spec  t  ivel  y  • 
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f'iR.  9 


SIMS  profile  of  300  keV,  U  x  10^^  cm*"^  ^^S  Implarkt  annealed  at  700*C 

for  20  mimit»*s;  Pearst)n  IV  profile  initial  eoiulition;  and  triple-zone 
profile  with  interfaces  at  0,A2  pm,  0.60  pm  and  I)|,  D2,  D3  equal 
^  K  lO'l^,  2.07  X  10“*^,  2,7  X  cm^/sec  respectively. 
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